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Abstract 

If X is a discrete abelian group and A a finite set, then a cellular automaton (CA) 
is a continuous map : J?' — >A^ that commutes with all X-shifts. If (/> : A — S'M, then, 
for any a € we define S(/>(a) = X^xex ^{'^yd (if finite); <p is conserved by ^ if Sc^ is 
constant under the action of J. 

We characterize such conservation laws in several ways, deriving both theoreti- 
cal consequences and practical tests, and provide a method for constructing all one- 
dimensional CA exhibiting a given conservation law. 

If ^ is a finite set (with discrete topology), and X an arbitrary indexing set, then (the 
space of all functions X ^) is compact and totally disconnected in the Tychonoff topology. 
If (X, +) is a discrete abelian groupQ (eg. X = Z^) with identity 0, then X acts on itself by 
translation; this induces a shift action of X on A^: if a = [oxlxgx] ^ -^^^ u G X, then 
<T"(a) = [6x|xex]> where = a(x+u)- 

A cellular automaton (CA) is a continuous map 5" : A^ — > which commutes with 
all shifts. The Curtis-Hedlund-Lyndon Theorem [Q] says that 3^ is a CA if and only if there 
is some finite B C X (a "neighbourhood of the identity") and a local map f : A^ — >A 

so that, for all a G A^ and x G X, 5^(a)x = f (ai ). Here, for any W C X, we define 
a|^ = [awlwGw] ^ ^'^^ example, if X = Z and B = [—B...B], then for any z G Z, 

^|b+z ~ ['^z-B, • • • , o-z+b]- Without loss of generality, we assume B is symmetric, in the sense 
that (b G B) (-b G B). 

^We assume X is abelian only for expositional simplicity; all results extend easily to nonabelian X. 
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Let +) be an abelian group (usually Z), and let : A — >Z. Heuristically speaking, 
(j){a) measures the "content" of a cell in state a E A. If is the identity of Z, we refer to 
= (;/)~"^{0} as the set of vacuum states, and 0^ = {a G A^ ; ax G 0, Vx G X} as the set 
of global vacuum configurations. If a G A^, then the function (f){a.) : X — >Z is defined: 
(/)x(a) = 0(ax). The support of a is the set supp [a] = {x G X ; ax ^ 0}; let A^^ be the 
set of elements of A^ with finite support. A CA ^ : A^ — > A^ is vacuum-preserving if 
;^(0^) C 0^, or, equivalently, if ^ (^<^) C A<^. 

Define E(/) : — >Z by: S(/)(a) = ^ 0x(a). li d ■ A'^ — > A^ is a CA, then we say 

xGX 

(f) is conserved by 5 if, for any a G A'^^, S(/)5^(a) = S(/)(a); we then write: (p G C{^; Z). 
Note that ^ must be vacuum-preserving to conserve (p. 

Example 1: Let A = Z = 'Ij/2, and let </> : A — >Z be the identity. If X = Z and 

oo 

a G A'^^, then S(/)(a) = an measures the parity of a as a sequence of binary digits. If 

n=— oo 

B = { — 1,0,1}, and f(a_i, ao, ai) = a_i + ao + ai, then 5^ is parity-preserving: T,(f)o^ = 
Ecj). 

Examples of Z- valued conservation laws for simple CA on (Z/2)^ are described in [Q. 
Necessary and sufficient conditions for conservation laws on one-dimensional CA are given 
in and used to completely enumerate the conservation laws for the 256 "elementary" 

(ie. nearest-neighbour) CA on (Z/2)^, and the 256 "elementary reversible" CA [|l^. 

Conservation laws arise most frequently in the context of particle-preserving cellular 
automata (PPCA). If 2 = Z and (p : A — then we interpret (j){ay.) as the number of 
"particles" at site x G X. Thus, S0 tallies the total number of particles in space; (p is 
conserved if particles are neither created nor destroyed. By extension, if Z = Z^ and 
(p : A — >N^ , then (p simultaneously tallies K distinct species of indestructible particles. In 
the simplest PPCA, cp : A — ^{0,1} (e.g. A = {0,1}, and (p is identity map); thus, at 
most one particle can occupy any site. PPCA on {0, 1}^ appear as models of traffic fiow 
|, |, |, [TT], 0, ig, and eutectic alloys |, ^. 

Example 2: If X = Z and B = { — 1,0, 1}, then there are exactly five PPCA with local 
maps f : {0, l}'^ — >{0, 1}. These are the identity map, the left- and right-shifts, and CA 
numbers 184 and 226 in the Wolfram nomenclature |T^. In CA #184, each "1" particle 



moves to the right whenever there is a "0" to its right, and remains stationary if there is 
a "1" to the right. CA #26 is the mirror image, with movement to the left. 

In §[l|, we characterize C{^; Z) in terms amenable to computational testing on a finite 
spatial domain. In §H, we use this to show how any (p : A — >M can be "recoded" by 
a M"'"-valued or N'^-valued function having equivalent conservation properties. In §& we 
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characterize C{^; M^) in terms of configurations with infinite support. In §^ we characterize 
CCS; R"*") in terms of spatial ergodic averages, assuming X is an amenable group, while in §^ 
we characterize C(^) in terms of stationary measures on A^, even when X is not amenable. 
In we consider the construction of CA with a particular conservation law. 



1 Finitary Characterizations 

Let O G be some fixed vacuum state. If V C X is finite and a G , then let (a) denote 
the configuration b G A^^ defined by: b|^ = a and = O for all x ^ V. Then define 

S0(a) = S(/)(a). If a G A^, then define S0(a)|^ = J] If U, V C X, define V + U = 

vev 

{v + u ; V G V, u G U}. If a G ^^+", then we define ^^(a) G by: ^^(a)^ = f (a|^^g) for all 
V G V. Let B(2) = b + 1; for example, if B = [-Bo...Bif C Z^, then B^^) = [-2Bo...2Bif . 

Proposition 3 G if and only if, for all a, c G A^^^\ identical everywhere except that 
tto = a c = Co, we have: 

5(c) - S(/) i?(a) = 0(c) -0(a). (1) 

Proof: Let 6 = 0(c) — 0(a). 

Proof of Clearly, S0((c)) = 5 + S0((a)), and thus, S05((c)) = 5 + S0i^((a)). 

Now, (a) and (c) only differ at 0, so ^((a)) |^^^ = d {{c}) |^^^, while i?((a)) |^ = ^ {a) 

and d{{c))\^ =d{c). Thus, 

S05(c)-S0i?(a) = S0i?((c))|^-S0^?((a))|^ = S0^?((c))-S0^?((a)) = 6, 

which yields equation (|I|). 

Proof of "<^=": Suppose a G A"^^, with supp [a] = {xi,X2, . . . ,X]\j}. For all n G [l..A^], 
suppose ax„ = and let 6n = 4>{hn)- 

Consider the vacuum a" G 0^ defined: 



O if x G supp [a] 
ay otherwise. 



We build a from a° one nonvacuum site at a time. For n G [l..A^], define a'^ G A'^^ by: 



hk if X = Xfc for some k < n, 
otherwise. 
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Thus, = a. For any n E [O..A^), a" and a"+^ differ only at x„, so 3^ (a") and ^J(a"+^) 
differ only in x„ + B; hence 



=[1] «) - <!> = ^n, 



n-l, 



where [1] follows from applying equation (0) at x„. Inductively, 5^ (a^) = 5„ + 

n=l 

E0 5(a°) = = S0[a] □ 



n=l 



Proposition ^ generalizes Proposition 2.3 of JT^ (which is the case X = Z). There is 
also a characterization of N-valued conservation laws in terms of periodic configurations (see 
Theorem 2.1 of |12] for case X = Z or Proposition 1 of |^ for case X = Z^); we generalize 
this to the following characterization of arbitrary conservation laws for any group X. 

Let be the set of all finite quotient groups X of X such that B*^^) maps bijectively onto 
its image B*^^^ C X under the quotient map. The local map f : — >A induces a cellular 

automaton ^ : — > ^^for any X G 0. Let C (^^^ = |0 : A — >M. ; (f) is conserved by . 
Corollary 4 <f) e C{d) iff (I) e C (d)- 

Proof: By Proposition H, G C{^) iff equation (P is true, while and G C (^^^ iff (|I|) 

is true, where (§ is (g) witj^"B(2)" replacing "B(2)". If we identify and in the 
obvious way, then clearly, (|ip and (Q) are equivalent. □ 

For example, if X = Z^, and B = [—B...B]^, then includes the quotient group 
X = CZ/m)^ for any M > A - B. Elements of correspond to M-periodic configurations in 
A^° (where a G J?''^ is M-periodic if <T^^'^(a) = a for any z G Z^). The action of ^ on J?- 
corresponds to the action of on M-periodic configurations in . Thus, ^ conserves iff 
conserves (j) on M-periodic configurations. 

2 Receding 

Proposition ^ yields a convenient "recoding" of real-valued conservation laws. Let = 
{x G M ; X > 0}, and let C {"S] K^) denote M+-valued elements of C (^^; M). 
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Proposition 5 Let cj) : A — s>M. 

1. There is a function (f> : A — so that, for any cellular automaton ^J, 

2. There is a function (f> : A — so that, for any cellular automaton ^, 

Proof: Part 1: Let — M = min^(a), and define <p(a) = (j)(a) + M for all a ^ A. Then 
clearly, satisfies the condition of Proposition |^ if and only if does. 

Part 2: A is finite, so (p^A) C M is finite, so the subgroup A C M generated by (p^A) 
is a finitely generated, torsion-free abelian group, therefore isomorphic to for some 
K. If C '■ ^ — is this isomorphism, and (p = C, o (f)^ then clearly {^(p E C{^] . 



We can always choose so that C l~l C N'^, and by Part 1 we can assume (p is 
nonnegative, so that (p : A — . □ 

Note that the vacuum states of (p are not the same as those of 0, because (p~^{Q} = 
(p~^{—M}. Thus, (p and (p determine two different notion of "finite support", and "^<^" 
refers to two different subsets of A^. 

Part 1 of Proposition 1^ implies that, to characterize real conservation laws, it is sufficient 
to characterize nonnegative ones; this will be useful in §^ and §^ Part 2 of Proposition ^ 
implies that we can interpret any real conserved quantity as tallying K species of indestruc- 
tible particles. Conversely, to construct a CA with a given real-valued conservation law, it 
is sufficient to construct a K-species PPCA; this will be useful in §0. 



3 A Nonfinitary Characterization 

Defining conservation laws in the context of A'^^ is somewhat unnatural, because A^^ is a 
very small subset of A^. We now characterize conservation laws in a way which is meaningful 
for any a G A^. 

For any W C X, define cl [W] = B + W, and int [W] = {w G W ; B + w C W}. Thus, 
W C int [cl [W]]. Because B is symmetric, int [W] = W \ cl [w"^] , where = X \ W. 
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Figure 1: 



Theorem 6 Suppose (p : A- 
any finite W C X 



". Then (j) G C(5') iff, for any a G with a' = 5^(a), and 



clP 



(2) 



Proof: Let V = cl [W] and let b = (^a|^^ —that is, b|^ = a|^, and 6x = o G 

for X ^ V, so b G A'^'^. Thus, if b' = S^(b), then b'|^ = a'|^, so it is sufficient to prove 
(Q) for b. The right-hand inequahty in (0) follows because: 

S<^(a')|w = S<^(b')|w <w S<^(b') S0(b) = S0(a)|^, 

where [1] is because is nonnegative, and [2] is because (j) G C(^). 

To see the left-hand inequality in (|), let Va = cl [V] and W = V2 \ W (see Figure |I|). 
Thus, supp [b'] C V2, so 



E0(b') 

If V3 = cl[V2], then cl [w 
E0(b) = 



m^X^ = S0(b')|^ + s0(b')|^, (3) 

V3 \ int [W], and clearly, supp [b] C W C V3, so that 



S0(b)| 



E<^(b) 



cl 



S0(b)| 



intP 



But applying the right-hand inequality in (H) to W, we have 



S0(b') 



< S0(b) 



cl 



while, by hypothesis that <p is conserved, we have 

E(/)(b) = S0(b') 

Combining yields: 



(4) 

(5) 
(6) 



cl 



+ S0(b) 



intr 



from which we conclude that S0(b) 



intr 



= + E^(b')|„ 

< S0(b')|^, as desired. 



"<^=": First, note that ^ must be vacuum-preserving: If a G 0^ and a' = ^(a), then (0) 
imphes that, for any x G X, < 0(aQ < S0(a)| =0, so a[, G 0. 

I cl [xj 

Next, suppose a G A<^, with supp [a] = Y. Let W = cl [Y]; then, since Y C int [W], (|) 
imphes: 

E0(a) = E0(a)|^ < < S0(a')|^ < S0(a)|^,j^, = E^(a), 

so that S^(a')|^ = S0(a). But ^ is vacuum-preserving, so supp [a'] C W; thus, 
S0(a')|^ = S0(a'), thus, S0(a') = E0(a) □ 



4 Conservation and Spatial Ergodic Averages 



A F0lner sequence [|T6[ on X is a sequence of finite subsets I„ C X so that, for any x G X 

card [In n (I„ + x)] 



hm 

n^oo card [I„J 



The group X is caUed amenable if it has a F0lner sequence. For example, is amenable. 



because I„ = forms a F0lner sequence. If X is amenable and Y C X, we define the 

Cesaro density of Y by 

, . , card[YnlJ 

density Y = lim , , ■ 7 

^ ^ J n^oo card [In] 

If a G then the (spatial) ergodic average of on a is defined: 

ErgAvex 0(a) = lim^y^0(ai) (where /„ = card [I„]). (8) 

ien„ 

If X = Z and I„ = [0...n], these correspond to the classical Cesaro density and ergodic 
average. We say that Y (respectively a) is stationary if the limit in (respectively, (H)) 
exists and is independent of the choice of F0lner sequence. 

Let A^°" [A^~\ be the set of probability measures on which are invariant under all X- 
shifts, and let M.'^ [^^] be the ergodic measures: the extremal points of A^°" [A^~\ within 
Ai [A^~\ . If /Lt G Ad'^ [^^] and X is amenable, then the generalized Birkhoff Ergodic Theorem 
says that yU-almost every a G A^ is stationary, and, if ;U G ^A'^ [^^] , then for yu-almost 
all a G A^, ErgAvCx 0(a) = (0, /i), where we use the notational convention: 



This yields the following characterization for conservation laws: 
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Proposition 7 Let X be an amenable group, : A — s>]R^. The following are equivalent: 

1- cm 

2. For any stationary a G J?', if sl = ^(a), then ErgAve^ 0(a) = ErgAve^ 0(a'). 

3. For any G M'^ [^^] , if fi' = ^{fi), then {(p, /i) = {(f), fi') . 

4. For any /i G [A^], if fi' = ^{fi), then (0, /i) = (0, ^'). 



Proof: (1^>2): Let be a F0lner sequence; if we define J„ = I^+M and K„ 

tlien {Sn}'^=i and {K„}^j^ are also F0lner sequences. Let /„ = card J„ = card [J„], 

and Kn = card []K„]. Since B is finite, the F0lner property implies: lim = 1 = lim — ^. 

n^oo J„ n^oo J„ 

Given e > 0, find ri G N so that 

/ K 

1-e < ^ and — ^ < 1 + e; (9) 

(1 - e) • ErgAvex 0(a) < ^$^0(ai); (10) 

I Gin 

^J]0(aO < (l + e)-ErgAvej,0(a). (11) 



Kn 

keK„ 



Now, In C int [J„] and K„ = cl [J„], so applying Theorem ^ to J„ yields: 

5^0(ai) < 5^0(a;) < 5^0(ak). (12) 

ien„ jejin kGK„ 

Thus, (1 - e)' ■ ErgAvex 0(a) <0 ^-i ^ 0(0,) <|| ^^(0,) 



\&n iei 



< 



jGJI,i 



< 



" keK„ " kGK„ 

<0 (l + e)2 ■ ErgAvex 0(a), 

where each inequality follows from formula with the same number. Letting e— s>0 as n^oo, 
we conclude by a squeezing argument: 

ErgAvex0(a) = lim(l±e)^-ErgAvex 0(a) = lim ^ ^ 0(an = ErgAve^ 0(a') 

jGjn 



(2=^1): Let b e ^<^, with b' = 5'(b). Let U = supp[b] and U' = cl[U], and find 

stationary Y C X with density [Y] = (5 > 0, such that yi + U' and y2 + U' are disjoint for 
any yi 7^ y2 G Y. Then define a G by: — «^^(b) for every y G Y, and Ox = O e 

for aUx^Y + U. Let a' = :J(a). ^ 

Claim 1: ErgAvex <^(a) = S ■ S(^(b) and ErgAvcx 0(a') = 6 ■ E0(b'). 

Proof: Let {lln}^i be a F0lner sequence; for any n e N, let Y„ = Y fl I„ and 
Y* = {y e Y ; y + U C Assume U contains 0, the identity element of X; thus 
Y* C Y„. By construction, for any y G Y, ^ (/» (cy+u) = E(/»(b). Thus, 

card [Y;] • E0(b) = J] Y.^M 

yeY* ueu 

< Xl-^^^i) ^ Yl Zl'^K+u) = card [Y„] • E(/)(b) 
iein yeY„ ueU 

Now divide everything by /„ = card and take the limit as n— >oo. By definition, 
density [Y] = 5, so that 

n— »oo n— >c» 

thus, by a squeezing argument, 

ErgAvex 0(a) = lim^y^0(ai) = ^•E0(b). 

" iei„ 

The proof for a' and b' uses U' instead of U, and the fact that, for any yeY, 

5]0(a;+,,)=^W ° [Claim 1] 

u'eu' 

Thus, S ■ E0(b) — ErgAvex 0(^) =[2] ErgAvex 0(^0 — ^ ' ^0(t>')) where inequality 
[2] follows from hypothesis (2). This implies E0(b) = E0(b'); since this holds for any 
b e A^^, we conclude that e C(g^). 

(2=^3): Apply the Birkhoff Ergodic Theorem. 

(3=^4): Any element of M."' [A^], is a weak*-limit of convex combinations of ergodic 
measures. So, suppose /x = wk*lim where Uj, — y with rjij e M.'^ [^^1 > 

i— »oo ' 

Xij e [0,1], for all i e N and j G If /u' = difJ-), then fi' = wk*lim z/-, where 
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y'i = dii^i) = ^>^ijVij, with rj'ij = diVij)- But by hypothesis (3), = {(i),Vij) for all 

i and j. Thus, (0, Ui) = (0, for all z G N; thus (0, /i) = {(/), /i'). 

(4^>2): If a G is stationary, let 5a G 7W [A^] be the point mass at a; then ^a' = 
di^a)- Let {I„}5^i be a F0lner sequence, and for all n G N, let /x„ = ^^^''"''^a and 

I6II7I 

fi'n = Y' ^^""'^a'- Since ^A [A^] is compact in the weak* topology, the sequence {fin}'^=i 



has a weak* cluster point, which by construction is shift-invariant. Dropping to a 
subsequence if necessary, we'll say /i = wk*lim Thus, 

n^oo 

= lim (0,/z„) = lim ^ (0, cr'^a) = lim^Y]0(ai) = ErgAve^ 0(a). 

iGln ieln 

If yu' = then fi' = wk*lim /i^ is also shift-invariant, and = ErgAve^^ 0(a'). 

ra— >oo 

But by hypothesis (4), we have (0,/^') = (0,/i); hence, ErgAve^ 0(a') = ErgAve^ 0(a). .□ 

Real-valued conservation laws thus preclude unique ergodicity: 

Corollary 8 Let Rq = min0(a) and Ri = max 0(a). //0 G C{^), then for any r G [Rq, Ri], 
there is an ^-invariant measure jJr G Ai'^ [^^] such that (0, fir) = r. 

Proof: Let G ^ be such that 0(afc) = Rk- Given r G [Ro,Ri], let A G [0, 1] be such 
that r = XRq + (1 — A)i?i. Let p be the probabihty measure on A with p{ao} = A and 
p{ai} = 1 — A, and let u^. be the associated Bernoulli measure on .4.^ — that is, the product 
measure = p. 

1 ^ 

Thus Ur is shift-ergodic and (0, Ur) = r. For all N gN, define rjN = — ^I'r', then by 

n=l 

part 3 of Proposition ^ (0, rjj^) = r. Since Ai'^ [^^] is compact in the weak* topology, 
the sequence {?^Ar}5^i has a weak* limit point, Pr- By construction, p^ is ^J-invariant, 
shift-invariant, and (0, /i^) = r. □ 
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5 Characterizations by Measure 



If X is not amenable, then the methods of §^ are inapphcable. However, we can still charac- 
terize conservation laws on in terms of shift-invariant measures, by projecting X onto a 
finite quotient group. 



Let n be as in 



If X G £}, then let M'' 



real- valued measures on A^, and Ai'^ 



A^ 



c M' 



A"^ 



be the space of shift-invariant, 
the space of shift-invariant 



probability measures. Let B C X be the (bijective) image of B under the quotient map 



X — >X, and, for any n E 



A' 



let 



prg*(/i) G M 



A^ 



be the marginal 



projection of /i onto B. Via the bijection 



I, we can identify A?" with AP", and thus 



M \A\ 



with M [A^; M] . Now define: 



M" [A^; M] = ^veM [A^; M] ; = /^^ for some /i G A^" J^] M and X G o}. 

Define ^jf) : A^__ — >M by 90 (a) = (f(a)) - 0(ao), for any a G J?. For any X G 0, define 
dt : A^ — by dt 0x (a) = (^^(a)^) - 0(ax), for all a G and x G X. 

Proposition 9 Let : ^ — anc? /ei ^ be a CA. The following are equivalent: 

1. 0gC(5). 

2. For allXe£}, (pec(d] 



3. For all X G £}, and all fx G M"" 



A' 



4. For all X G 0, and all fi G M"" 

5. For all fig G M" [A^; M] , (50, fig) = 0. 
Proof: (HI) <^ (H): This just restates Corollary | 



, (^4 00, fl) = 0. 



^ (D: If fieM'' 



A' 



, then (9*00, /i) = (50, /ig). 
i) ^ (i: By definition, {dt 0o, /i) = (0o o - 0o, /i) = (0o, d{f^)) - (0o, m)- 
(D^d): X is finite, so = A^. Thus, we can well-define : 



by: 



9i0(a) = ^at0x(a). If ^ G 

(5i 0,/^)- 



and we likewise define fi = cr^/^, then 



xGX 

clearly, (5f 0o, /i) 
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But if (/) G C ( , then dt()) = 0. To see this, let a G A^; then: 
= 5^(0x^?(a)-</)x(a 



S0;^(a)-S0(a) 



0. 



xex 



Also, if /i G A , Rj , then = card X ■ /i. Combining these facts yields: 
card X -(9*00, A*) = (5* 0o, Ji) = (dtcp^fi) = (0,/^) = 



0. 



(1^): If a G and <5a G A^ is the point mass at a, then 6^ = cr'^'^a is in 



. Thus, 

E0^?(a)-E(/*(a) = (0, 5^(a) - 0,(a)) 

xex 

= (5a) = 

This is true for any a G A^, so (p E C [^) . 



{dt(J),T^) = 0. 



□ 



A^ is finite, so Ai [A^; M] is a finite dimensional vector space, and [A^; M] is a 
linear subspace, with some finite basis Hi, . . . , fiN (for example, see [0). To check Part 5 
of Proposition ^ it suffices to check that (90, = for all n G [l..A^], a finite system of 
linear equations. 

Corollary 10 Let card [A] = A and card [B] = B. If cj) e C(S^), then ^ o f(a) = 
A^^^''*^(j){a). In particular, if A = [0...A) C N, and (j){a) = a for all a E A, then 



aeA 



-A^{A-1). 



Proof: Let r] G Ai'^ [A^] be the uniform Bernoulli measure, assigning probability to 
every element of A^. Then o f(a) — - 0(a) = (0 o f , r^) — (0, 77) = 



(90, r/) 



0, by Part 5 of Proposition ^. Thus, o f (a) = — 



fa). _□ 



The second statement of Corollary |T0| generalizes Corollary 2.1 of which is the 
analogous result when X = Z and B is an interval. 



12 



6 Constructing Conservative Cellular Automata 



Given (j) : A — can we construct a cellular automaton 3" : — ^ so that G C{^)1 
Can we enumerate all such automata? By Part 2 of Proposition we can assume : 
A — ^'H^; hence, the problem is to construct a CA which preserves K species of particles. 
For simplicity, we will consider the construction of a CA preserving one species of particle. 

PPCA are usually constructed by explicitly specifying how each particle displaces from 
its current position to a nearby location, in response to its current local environment. The 
local rule of the PPCA, as a map f : — >A, is then formulated a posteori to realize this 
"displacement" model. Does every PPCA arise in this manner? Or are there PPCA not 
admitting any displacement representation? 

In the case X = Z, ^ = [0...74], and 0(a) = a, every PPCA unearthed in extensive 
computational searches has admitted a displacement representation |12[. But it is not 
clear why such a representation should always exist, and, for more complicated PPCA, such 
a representation, even if it exists, may not be obvious from inspection. 

If a displacement rule is to yield a cellular automaton, it must satisfy the following 
conditions: 

(Dl) The rule is equivariant under shifts: a particle at x G X, in configuration a G A^ will 
experience the same displacement as a particle at x + y in configuration cr^(a). 

(D2) Each particle has bounded velocity. The new position of any particle at x is inside 
x + B. 

(D3) Each particle's displacement locally determined: The displacement of any particle 
at X is entirely determined by a| ^^^^^ . (Heuristically speaking, when "deciding" its 
trajectory, a particle must look not only in a B-neighbourhood around its current 
location, but also in a B-neighbourhood around each of its possible destinations.) 

Several particles may be present at a given site; the rule must assign a displacement to 
each of them, yielding a multiset of displacements, which can be represented as an element 
of N^. Formally, a particle displacement rule^ (PDR) is a function : A^^''^ — m*. If 

a G A^, then for all x G X, we write Ox(a) = ^a| ^^^^^ j and write the components of this 

object as Ox^y(a) for all y G (B + x). 

If G C{^; N), then we say that is compatible with and ^ if, for any a G A^^, with 
a' = ^^(a) the following two conditions hold: 

(CI) 0(ao) = 5^c)o^y(a) and (C2) 0(0 = 5^0x^o(a). 

yeB xeB 

Given and d, we can construct all 0-conserving CA compatible with d as follows: 
^This generalizes the particle automata introduced in lit]. 
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1. Let M = max(^(a). For every m G [0..M], let Am = 4>^^{m} C A. 

2. For any a G A^^''^\ if m = c)x^o(a) then define f(a) to be any element of Am- 

Thus, if card [Am] < 1 for all m, then the PDR uniquely determines f, but if card [Am] > 
1 for some m, then there are several CA compatible with d. Heuristically speaking, if 
card [Am] < 1 for all m, then particles have no 'internal state', and are totally interchange- 
able. If card [^i] > 1, however, then a solitary particle has more than one 'internal state'; 
if card [Am] > 1, then a pile of m particles at a site can be arranged in more than one 
'configuration'. The PDR assigns locations to the particles, but the local rule f is still 
responsible for deciding what 'configurations' and 'internal states' they assume. 

Note that, in general, this method yields a local rule f : A^''^^ — >A. However, if d is 
compatible with a cellular automaton ^ having local rule f : A^ — >A, then condition (C2) 

ensures that c)x^o(a) will also depend only on aj^^ ; thus, f will be a function f : A^ — >A. 

The problem of constructing a PPCA is thus reduced to the problem of constructing 
a PDR. To show that every PPCA arises in this manner, it suffices to show that every 
PPCA has a compatible PDR. For arbitrary X, this is surprisingly difficult; there are many 
directions a particle can go in, and potentially several particles vying for each destination. 
When X = Z, the one-dimensional topology obviates these complications. 

Proposition 11 If ^ : A^ — > A^ is a cellular automaton, and (j) G C(^;N), then there is a 
PDR compatible with ^ and □ 

We will construct the PDR via a naturally defined "flux" function, which describes the 
flow rate of particles past each point in Z. Assume B = [—B...B]. If z G Z, then we deflne 
the flux from z to z + 1 as follows. Let a G A'^'^ and let a' = 5^(a). Since is conserved, we 
know that 

z oo oo oo z oo 

y=— oo y=z+l y=— oo y=— oo y=— oo y=z+l 



Thus, 



^ <^(ay) - 0(4) = ^ 0(4) - ^ 0(ay) (13) 

y=— oo y=— oo y=z+l y=z+l 



Let Iz (a) be the quantity on either side of (p!3|); this is the flow from z to z + 1. Thus, 

Iz (a) = — Iz-i (a) is the flow from z to z — 1, and Iz (a) = Iz (a) + Iz (a) is the total flux 
out of site z. 
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Proposition 12 Let a e and a' = ^^(a). //0 e C{^; N), then: 

1. Iz (a) = -dtcpz (a). 

2. For any z, the value of Iz (a) is a function only o/aj^^^^ . 

z z+B 

3. (i) Iz(a) < '^(«y)' I^(^) ^ SZ-^^)' 

y=z-B y=z 
^ z+B ^ z 

(^^^) Iz (a) < J^^K); ^ Iz (a) < Yl ^^^4) 

y=z y=z—B 

Proof: Without loss of generality, assume z = 0. 
Part 1 follows from equation (|T3|) by straightforward algebra. 



Part 2: Suppose b G with aL „ = bL „ . Define c G by cL = 

|[-_B..._B] |[-_B..._B] \{—oD...B] 

aL while cL ^ = bL „ ,. Thus, if b' = ^?(b) and c' = 5(c), then c'l, „ = 

I ( — CO.. .is] |[— i!...oo) I [— i;...oo) ^\ / ^\ // |( — OO...UJ 

^'|(-oo...O] ^^il^ ^'|[0...oo) ='^'|[0...oo)- Thus, 

Io(a) = ^ <^(ay) - ^ <^(a;) = ^ 0(cy) - ^ 0(4) = Iq (c) 

y<0 y<0 y<0 y<0 

0<y 0<y 0<y 0<y 

Part 3: To prove (i) let b = (ai \, Thus, Iq (b) = Iq (a) by Part 2, but 6y = o G for 



y 



all y ^ B, so that lo (b) < ^ '^f^^) = E '^(^y) = E "^("y)' 

y<0 y=-B y=~B 

Inequality (iii) is proved similarly. For inequality (ii), suppose that supp [a] U supp [a'] C 
[-A^...oo]. Then 



lo (a) = E '^(%) " E '^'^4) = E '^'^"y) " E '^'^"y) 

y<0 y<0 y=-N y=-N-B 

OB B 

= E '^(«y)- E '^K) + E'^K) 

y=-N y=-N-B y=l 

B B 

<[i] Y (t^i^y)- E '^(^y) + E'^^^) = E'^K)' 

y=-Af y=-N y=l y=l 

where [1] follows from Theorem^, with W = [— A^...O]. Inequality (iv) is similar. □ 
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Part 2 of Proposition |T^ implies that we can well-define Iz (a) even when a has infinite 
support. Note that Part 2 does not imply that 

B B 

Io(a)= 0(ay)- '^K) = Xl'^K) ~ Xl'^K)- 

y=-B y=-B y=l y=l 

Example 13: 

(a) Suppose = <tMs the five-fold shift. If a = [. . . 10111100110100101 . . .], then a = 
[...10111100110100101...], and Iq (a) = 3 (here, the hat indicates the position of 
ao). If a = [. . . lllllilllll . . .] , then Iq (a) = 5. 

(b) Suppose ^ is CA #187 from Example |[ If a = [. . . 10 . . .] , then Iq (a) = 1, while if 
a = [. . . 11 . . .] , then Iq (a) = 0. 



Proof of Proposition 11: Again, assume without loss of generality that z = 0. We 



construct the PDR from the flux as follows: 

Case 0: (Iq (a) < and Iq (a) < 0) No particles leave site in either direction, so 
fo(a) = 0. 

Case 1: (0 < Iq (a) and < Iq (a)) Now, particles are leaving site in both directions. 
For z > 0, if (f>{o,y) < Iq (a) then define 

0<y<z 



c)o^z(a) = min 



Io(a)- J2 '/'K)[ 

I 0<y<z J 



and set 5o^z(a) = if Iq (a) < ^ 0(4)- Likewise, for z < 0, if ^ 0(4) < lo (^) 

0<y<z z<y<0 



then define 



c)o-.z(a) = min <^ 0(4), lo (a) - ^ 



l«yJ 



z<y<0 



and set c)o^z(a) = if lo (a) < ^ ( 

z<y<0 



ia': 



Case 2.1: (lo (a) < < lo (a)) Now, particles enter from the left, and leave to the 
right. Let Jo = 0(ao)) fi'^st Jo particles entering from the left will fill the Jo available 
destinations at 0. Let Ji = max{0, — lo (a) — Jo}; the next Ji particles entering from 
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the left must pass through 0, and will fill up the next Ji destinations available to the right 
of 0. 

Iq (a) particles leave to the right. If Ji > 0, then the first Ji of these particles are from 
(— OO...0), while the last Iq (a) — Ji are the particles originating at itself (Part 1 of 

Proposition |l^ implies that Iq (a) — Ji = Iq (a) + I o (a) + (p^tt'o) = —dt 4>o (a) + (;/>(aQ) = 
0(ao)). 

Thus, for all z > 0, if Ji < < lo (a) then define 

0<y<z 



^o^zi^) = min <^ lo (a) - ^ 0(a^ 



0<y<z 

and set ()o^z(a) = if 0(<3y) < Ji or if Iq (a) < 0(fly)- This leaves the boundary 

0<y<z 0<y<z 

case when 0(Oy) < Ji < 0(fly)- In this case, let c)o->z(a) = (pWy) ~ Ji- 

0<y<z 0<y<z 0<y<z 

Case 2.2: (Iq (a) < < lo (a)) Now, particles enter from the right, and leave to the 
left. This case is handled analogously to Case 2.1 

It remains to verify that c) is a PDR, and is compatible with 5^. In each of Cases 1, 2.1, 
and 2.2, Do-+z(ci) = for any z ^ B; this follows from inequalities (ii) and (iv) in Part 3 of 
Proposition |l^. Thus, condition (D2) is satisfied. Also, the value of fo(a) is determined by 

a'|j^, Iq (a) and Iq (a). Clearly, a'j^^ is determined by aj^^^) i while Part 2 of Proposition 

12l says that Iq (a) and Iq (a) are determined by aj^^. Thus, c)o(a) is a function only of 
ai , so condition (D3) holds. We have described the algorithm at z = 0, but we apply 
the same algorithm at all points, so condition (Dl) holds automatically. 

Compatibility condition (CI) holds by construction in each Case. To check (C2) suppose 
a G v4.^^, with supp [a] = [— A^...A^]. Beginning at — and proceeding to the right, we can 
inductively verify (C2) for each z G Z by applying Part 1 of Proposition □ 



7 Conclusion 

In §|I] and we provide practical methods for detecting the existence of conservation laws, 
while in §^ and we provide abstract characterizations of such laws. Theorem ^ from §^ 
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is used in §| to develop a notion of 'flux', which yields a method for constructing cellular 
automata having a particular conservation law. 

However, the 'displacement representation' constructed in §^ is inapplicable to the case 
X = iP ^ D > 2. Do displacement representations exist for PPCA on higher dimensional 
lattices? If we interpret a conserved quantity as the density of some material, many ques- 
tions remain about the 'hydrodynamics' of this material: its patterns of flow, concentration, 
and diffusion. We also expect that higher-dimensional PPCA may exhibit complex parti- 
cle dynamics, including the formation of complex, large-scale, stable clusters analogous to 
molecules. What is a good framework for studying these quasichemical dynamics? 
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